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We present a spectral analysis of the central X-ray emission for a sample of galaxy 
clusters observed with Chandra. We constrain the quantity of a second cospatial 
temperature component using Markov Chain Monte Carlo sampling and discuss the 
implications for our understanding of cooling flows. 



1. Introduction 

The cores of many galaxy clusters are sufficiently dense 
and cool that the plasma cooling time is shorter than a 
Hubble time. For many years it was thought that run- 
away coo ling would result in a large central mass deposi- 
tion rate llFabiarJl994HAllen FabiarJl997tlPeres et all 
Il998t IWhitdl2nnnHAllenll20n(TD . Chandra and XMM ob- 
servations have altered this picture significantly - there 
appears to be a core temperature floor of 1-2 keV, and 
inferred mass depos ition rates have been reduc ed by an 
order of magnitude l|Peterson et alJl200ll [2003) . Several 
mechanisms have been proposed to explain the lack of 
colder gas, including heating by AGN, heat conduction 
from cluster halo plasma, and small-scale variations in 
the cooling and metallicity structure of the plasma. Each 
of these processes can potentially leave a specific obser- 
vational signature. For example, if conduction provides 
the energy required to arrest the cooling then the core 
plasma may be single-phase, but if AGN and/or small- 
scale inhomogeneities are responsible, one might expect 
to see observational signatures indicating the presence of 
multiphase plasma. 

We have analyzed a sample of 12 galaxy clusters found 
in the data archive of the Chandra X-ray Observatory for 
evidence of multiphase plasma in each cluster core. We 
first describe the method, briefly sketching the Markov 
Chain Monte Carlo technique, and apply it to a sample 
of clusters. We then discuss the results and their impli- 
cations for cooling flows. 
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Fig. 1. — Core-halo geometry of the simple and complex models. 
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Fig. 2. — Schematic representation of the simple and complex 
model parameter spaces. 



2. Method 

We use a simple core-halo geometry and Markov Chain 
Monte Carlo (MCMC) simulations to assess the statis- 
tical significance of the presence of multiphase plasma 
in a sample of clusters. We compare two models of the 
emssion: a simple model M s which contains one emission 
component in the halo and one in the core, and a com- 
plex model M c which contains one in the halo and two 
in the core (see Figure 

Each model is fit to a data set consisting of two spec- 
tra, one from the outer annulus and one from the inner 
annulus. Because M s lies on a boundary of the parame- 
ter space of M c (i.e. the normalization of the second core 
component goes to zero; Figure EL th e standard F test 
cannot be used l|Protassov et al.ll2002|) . We must instead 
construct an empircal F distribution, using MCMC sam- 



pling (sketched below), to which we apply the F test. 
The F statistic, 

x 2 (M s |D)-x 2 (M c |D) 

X 2 (M S |D)/KM S ) ' 1 ' 

is a measure of the improvement in a fit when the simple 
model is replaced by the complex model. The F distri- 
bution of a large number of data realizations allows us to 
quantitatively compare these two models by comparing 
the F distribution to the F value of the original data set. 

Given the (known) probability distribution function 
P(x), where x is the vector of model parameters, we 
construct an empirical F distribution using MCMC sam- 
pling and perform an F test to determine the significance 
of a second emission component in the core. The entire 
procedure is as follows: 
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1. Model the real data set Do with M 5 using 
XSPEC; call the best-fit parameters Xg. 

2. Use XSPEC to calculate P(D |x) (i.e., the 

likelihood) . 

3. Use Bayes' Theorem to calculate P(x|D ). 

4. Create a large sample of model parameters 
x^ using a random walk through through 
the parameter space, rejecting each new 
position which does not meet an accep- 
tance criterion based upon P(x|D ). This 
is the Metropolis algo rithm form of the 
MCMC technique; see lNeall (119931) . (We 
also discard all unphysical excursions in 
parameter space, i.e. where T < or 
p < 0.) For each x|, compute a fake data 
set Di, including instrumental effects of 
the Chandra telescope and detectors, as 
well as counting statistics. 

5. Fit each M s and M c to each D 4 . 

6. For each pair of models tabulate its F 
value given by Equation ^ 

7. Bin up the set of F values, creating an 
unnormalizcd histogram , and compare to 
the F value of the original data. 

In many applications of MCMC sampling one pays spe- 
cial attention to the finite "burn-in" period during which 
the Markov chain equilibrates. The length of the burn-in 
phase depends upon the sensibility of the starting point, 
and the appropriateness of the scale chosen for the pro- 
posal probability distribution step. This is not a consid- 
eration in our case because we start each MCMC sample 
at the (already known) peak of the probability distribu- 
tion function -P(x). 

In practice this recipe is computationally intensive, not 
because of any features of the MCMC sampling per se, 
but because each of the faked spectra must be modelled 
twice. For a sample size of 1000 simulations, XSPEC 
must simulate 1000 spectra and calculate 2002 sets of 
best-fit values for the model parameters (including the 
original data). This fact leads us to simplify the method. 
First, in order to reduce the modelling time, we have 
adopted a simplified core-halo geometry. In this scheme 
the "core" is represented by a single shell (in this case a 
sphere), while the halo is represented by another shell. 
Thus M 5 contains four parameters, the temperature and 
density of each of the two shells, while M c contains six, 
the additional two parameters representing the temper- 
ature and density of a second cospatial emission compo- 
nent in the core. This simplification also greatly improves 
the numerical stability of the fitting procedure. The al- 
gorithm (steps 1-6 above) is implemented in a Tel script 
run within XSPEC. 

Once we have completed step 7 we can distinguish be- 
tween the models. The location of the F value of the 
data within this empirical F distribution contains infor- 
mation regarding the relative merit of M s and M c . We 
define the significance S of the distribution as 



J ^N(F)dF 
S J oo N(F) dF [Z) 

The signficance S = 1 — Pf , where Pf is the probability 
that the simple model constitutes the better description, 
and that the F value of the data is this large strictly by 
chance. Thus, for a one-parameter model, S — 0.68, 0.90, 
and 0.99 may be interpreted as 1, 2 and 3tr detections of 
the additional component. When discussing the presence 
of a second emission component in the X-ray data, we 
adopt 99% as a threshold significance. 

3. Application to Chandra Spectra and Clusters 

We apply this method to 12 clusters observed with 
Chandra. All clusters in the sample are fairly round and, 
with the exception of CL0024, seem to be fairly spherical 
with no significant amounts of substructure. Ten of these 
clusters are generally agreed to contain "cooling flows" in 
their center s , whil e two of them, CL0024 and A2104, do 
not. (200QD derives an upper limit to the cooling 

rate in A2104.) In each cluster, the core size is defined 
as the radial extent containing roughly 3000-4000 source 
photons. Except for Hydra A, all of the clusters in the 
sample admit a second emission component in this re- 
gion. (In Hydra A, the best-fit temperature of the second 
component equals that of the first, rendering an MCMC 
F-test irrelevant.) 

The resulting F distributions for the entire sample are 
shown in Figure Statistics for each cluster are shown 
below. We list list the cluster, the pre-Chandra/XMM 
cooling rates from the literature, a reference for this 
value, and the multiphase plasma MCMC signficance. 

Table 1.— The 12 Chandra clusters in the sample. 
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cluster 


(M y 1 ) 


reference 


(N = 10C 


A1689 


118 ±f Ig 


White 


(2000) 


34.2% 


A1795 


453 ±H 


White 


(2000) 


82.3% 


A1835 


683 ±f 7 l ! 


White 


(2000) 


48.3% 


A1942 


817 ±£if 


White 


(2000) 


17.6% 


A2029 


547 ±l\ 


Allen 


(2000) 


99.8% 


A2104 


Q±f 


White 


(2000) 


66.1% 


A2204 


984 ±1|| 


White 


(2000) 


99.9% 


CL0024* 








77.0% 


HydraA 


264 ±%l 


Allen 


(2000) 




MS1358 


691 ±1|? 


Allen 


(2000) 


63.4% 


MS2137 


1467 ±?|g 


Allen 


(2000) 


27.1% 


ZW3146 


2228 ±@36 


Allen 


(2000) 


98.9% 



* There are no measurements (nor upper limits) of 
cooling flow plasma in CL0024 in the literature. 



4. Discussion 

Of the 10 cooling flow clusters in the sample, only three 
of them - A2029, A2204, and ZW3146 - show evidence 
for multiphase gas in the core, if we adopt S = 99% 
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Fig. 3. — Empirical F distributions for the 12 clusters in this 
study. Note that Hydra A does not admit a second core component 
whose temperature differs from the first. 



as our significance threshold. Perhaps more surprising 
is the fact that the two clusters with the largest pre- 
Chandra/XMM mass deposition rates - MS2f37 and 
ZW3146 - display such dramatically different evidence 
for the existence of multiphase plasma. Both were sus- 
pected of harboring cooling flows with rates in excess 
of fOOO M Q y" 1 , and yet MS2137 shows no evidence of 
multiphase gas. (As expected, the two non-cooling flow 
clusters also show no evidence for a second core emission 
component.) ft could be that, in the absence of recent 
core merging events, the equilibrium state of the plasma 
is uniphase, and that the merger and accretion of smaller 
(and cooler) substructures during the continuing assem- 
bly of the cluster are responsible for some clusters show- 
ing evidence of multiphase cores. If this were the case, 
one might look for some evidence for a merging event in 
MS2137 which has occured within a cooling time. (As 
yet we have not found none.) Regardless, it seems that 
cooling flow cluster cores are an inhomogeneous class. 
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